Search for: All records

Abstract This study integrates data from all broadband seismic stations in Alaska and northwestern Canada in 1999–2022 to construct a shear‐wave velocity model for south‐central Alaska and northwesternmost Canada, using ambient noise wave propagation simulation and inversion. Our model reveals three key features, including (a) the presence of the subducting Yakutat slab with apparent velocity reductions near the trench and within its flat segment, (b) two slab segments beneath the Wrangell volcanic field, differing in steepness, depth, and seismic velocity, and aligning spatially with the northwestern and southeastern volcano clusters, and (c) the existence of slab windows between the Yakutat and Wrangell slabs and between the northwestern and southeastern portions of the Wrangell slab. Our findings reinforce that the Wrangell volcanoes are predominantly influenced by subduction‐related magmatism. Furthermore, the two slab windows could have induced asthenospheric upwelling, contributing to the volcanism in the Wrangell clustered volcanoes. 

A typical subduction of an oceanic plate beneath a continent is expected to be accompanied by arc volcanoes along the convergent margin. However, subduction of the Cocos plate at the Middle American subduction system has resulted in an uneven distribution of magmatism/volcanism along strike. Here we construct a new three-dimensional shear-wave velocity model of the entire Middle American subduction system, using full-wave ambient noise tomography. Our model reveals significant variations of the oceanic plates along strike and down dip, in correspondence with either weakened or broken slabs after subduction. The northern and southern segments of the Cocos plate, including the Mexican flat slab subduction, are well imaged as high-velocity features, where a low-velocity mantle wedge exists and demonstrate a strong correlation with the arc volcanoes. Subduction of the central Cocos plate encounters a thick high-velocity feature beneath North America, which hinders the formation of a typical low-velocity mantle wedge and arc volcanoes. We suggest that the presence of slab tearing at both edges of the Mexican flat slab has been modifying the mantle flows, resulting in the unusual arc volcanism. 

The Cascadia subduction zone, where the young and thin oceanic Juan de Fuca plate sinks beneath western North America, represents a thermally hot endmember of global subduction systems. Cascadia exhibits complex and three-dimensional heterogeneities including variable coupling between the overriding and downgoing plates, the amount of water carried within and released by the oceanic plate, flow patterns within the mantle wedge and backarc, and the continuity and depth extent of the subducting slab. While recent research has benefitted from extensive onshore and offshore deployments of geophysical instrumentation, a consensus on many important aspects of Cascadia’s magmatic, tectonic, and geodynamic setting remains elusive. 

Abstract The northwestern part of North America has recorded multiple tectonic events, such as terrane accretion, strike‐slip motion, and subduction of the Pacific and Yakutat plates, providing an iconic setting to investigate the tectonic evolution of the continental crust. In this study we analyze the receiver functions at seismic stations deployed during 1999–2022 to estimate the crustal thickness, as well as possible slab signature, in Alaska and northwestern Canada. The Moho signal can be clearly detected within the continental region. Specifically, in northwestern Canada, the thickest crust is observed beneath the Cordilleran Deformation Front, which marks the structural boundary between the North American Craton and the North American Margin. We observe a few distinct offsets in the Moho depth located both within the tectonic units and approximately across the major faults between the tectonic units. We provide a first‐order estimate of the depth gradient of the Moho offsets based on the horizontal distance of the two closest seismic stations across the offsets. We propose that the Moho offsets reflect the cumulative impact of the accretionary orogenies and post‐orogenic tectonic events on crustal modification. The continental Moho signal is weak or obscure in Aleutian and southcentral Alaska, and the oceanic Moho within the subducting plates is likely detected. This study provides new seismic insights into understanding the impacts of the tectonic events on continental formation and evolution. 

Abstract Subduction of the Nazca plate results in the uneven distributions of earthquakes and arc volcanoes along the South America's western margin. Here, we construct a high‐resolution shear‐wave velocity model from immediately offshore to the backarc in South America, using advanced full‐wave ambient noise tomography. Our new model confirms and provides further constraints on three major features, including (a) extensive low‐velocity anomalies within the continental crust, (b) two high‐velocity flat slab segments located beneath southern Peru and central Chile, and (c) complex slab geometry at flat‐to‐normal transitional subduction. The flat slab segments roughly correlate with the volcanic gaps but not with the seismicity gaps. We suggest that variations of slab geometry along strike and down dip have significantly modified the flow patterns within the mantle wedge. Subduction of oceanic ridges has altered the slab dehydration processes, which can influence the distribution of arc volcanism and the occurrence of intermediate‐depth earthquakes. 

Abstract The composition of the lower continental crust, as well as its formation, growth, and evolution, remains a fundamental subject to be understood. In this study, we carry out a comparative and integrative analysis of seismic tomographic models, teleseismic receiver function results, and Airy isostasy in order to investigate the properties of the lower continental crust in eastern North America. We extract the depths for Vs of 4.0 km/s, 4.2 km/s, and 4.5 km/s from three selected tomographic models and calculate the differences between the Vs depth contours and the Moho depth defined by receiver functions. We then calculate the Airy isostatic Moho depth and its misfit with the receiver‐function‐defined Moho. Our analysis reveals three key features: (a) the deepening of the Vs depth contours and the strong negative Airy misfit within the U.S. Grenville Province; (b) a seismically faster‐than‐average and compositionally denser‐than‐average lowermost crust in the eastern North American Craton and the Grenville Province; and (c) the thickest, seismically fastest, and densest lowermost crust beneath the southern Grenville Front, the southern Grenville‐Appalachian boundary, and the U.S.‐Canada national border. We suggest that the lower crust of the craton and the Grenville Province has densified through garnet‐forming metamorphic reactions during and after orogenesis, contributing to the widely distributed fast‐velocity layer. The lower crust beneath the tectonic boundaries could have experienced more extensive garnet growth during orogenesis and emplacement of mafic magma. This study provides new constraints on the seismic and compositional properties of the lower crust in eastern North America. 

Significant along-strike variations of seismicity are observed at subduction zones, which are strongly influenced by physical properties of the plate interface and rheology of the crust and mantle lithosphere. However, the role of the oceanic side of the plate boundary on seismicity is poorly understood due to the lack of offshore instrumentations. Here tomographic results of the Cascadia subduction system, resolved with full-wave ambient noise simulation and inversion by integrating dense offshore and onshore seismic datasets, show significant variations of the oceanic lithosphere along strike and down dip from spreading centers to subduction. In central Cascadia, where seismicity is sparse, the slab is imaged as a large-scale low-velocity feature near the trench, which is attributed to a highly hydrated and strained oceanic lithosphere underlain by a layer of melts or fluids. The strong correlation suggests that the properties of the incoming oceanic plate play a significant role on seismicity. 

Abstract The eastern North American passive margin was modified by Mesozoic rifting. Seismic data from recent deployment of onshore and offshore stations offer a unique opportunity for studying the signature of syn‐rifting and postrifting in lithospheric structures. Using full‐wave ambient noise tomography, we construct a new seismic velocity model for the lithosphere of the southeastern United States. Our model confirms an oceanic‐continental transitional crust over a ∼70 km wide zone across the coastline. Our model reveals (a) a patch of lower‐than‐average mantle lithospheric velocities underlying this transitional crust and (b) a low‐velocity column in the mantle lithosphere beneath the Virginia volcanoes. We propose that anomaly 1 represents cooled enriched mantle that underplated the thinning crust during the initial stages of rifting around 230 Ma. Anomaly 2 likely has a more recent origin in the Eocene and may result from an asthenospheric upwelling induced by a localized lithospheric delamination. 

Abstract Lithospheric layering contains critical information about continental formation and evolution. However, discrepancies on the depth distributions of lithospheric layers have significantly limited our understanding of possible tectonic connections among the layers. Here, we construct a high‐resolution shear velocity model of eastern North America using full‐wave ambient noise simulation and inversion by integrating onshore and offshore seismic datasets. Our new model reveals large lateral variations of lithosphere thickness approximately across the major tectonic boundaries, strong low‐velocity anomalies underlying the thinner lithosphere, and multiple low‐velocity layers within the continental lithosphere. We suggest that the present mantle lithosphere beneath eastern North America was formed and modified through multiple stages of tectonic processes, among which metasomatism may have significantly contributed to the observed intralithospheric low‐velocity layers. The sharp thickness variation of lithosphere promoted edge‐driven mantle convection, which has been consequently modifying the overlying mantle lithosphere and further sharpening the gradient of lithosphere thickness